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ABSTRACT 

The  U.S.  Navy  is  currently  studying  the  use  of 
double  hull  designs  in  high  strength  low  alloy 
(HSLA)  steels  for  surface  combatant  ships.  A  full 
scale  prototype  double  hull  module  was  fabricated, 
from  which  multicellular  box  column  specimens 
were  cut  for  compressive  tests  to  failure.  Initial 
imperfections,  i.e. ,  initial  plate  deflections  and 
welding  residual  stresses,  affect  the  stiffness  and 
strength  of  welded  members.  This  paper  describes 
the  measurement  of  these  imperfections  and  the 
analysis  of  their  effects  on  the  component  plates  of 
the  cellular  box  specimens.  Initial  deflections  were 
measured  in  the  laboratory,  where  the  maximum 
values  did  not  exceed  the  Navy's  guidelines  or 
proposed  values  of  several  researchers.  Residual 
stresses  in  a  box  specimen  were  also  measured  in 
the  laboratory  under  more  controlled  conditions. 
Using  the  measured  imperfections,  the  plate 
arrangements  were  analyzed  using  the  finite 
element  method.  The  imperfections  were  found  to 
reduce  the  stiffness  and  strength  of  the  plates.  The 
results  showed  that  for  accurate  prediction  of  the 
strength  of  welded  plates,  initial  imperfections 
must  be  taken  into  account. 

INTRODUCTION 

Many  double  hull  ship  designs  consist  of  twin 
skins  of  plating  wrapping  around  the  ship  bottom, 
sides  and,  optionally,  the  main  deck.  The  nearly 
parallel  plates  are  separated  and  stiffened  by 
longhudinal  web  girders  that  span  between 
transverse  bulkheads.  Other  transverse  components 
are  eliminated,  which  creates  a  simple 
unidirectional  structure  in  the  longitudinal 
direction.  The  U.S.  Navy  is  currently  studying  the 
implementation  of  an  advanced  double  hull  design 
for  surface  combatant  ships  (1) .  In  the  commercial 


sector,  product  oil  earners  and  tankers  have  already 
been  built  in  Japan  and  Korea  (2,3).  Some  of  the 
reasons  for  considering  double  hull  ship 
construction  include:  a  simplified  and 
unidirectional  structural  system  that  facilitates 
automated  welding  and  fabrication  technigues; 
inherent  strength  for  improved  ship  survivalibility 
in  collision  or  combay;  and  fewer  areas  of 
discontinuities  and  complex  welded  details  that 
could  give  rise  to  fatigue  and  fracture  problems. 
Figure  1  shows  the  difference  between 
conventional  and  double  hulls  of  ships. 


Figure  1:  Conventional  and  advanced  double 
hull  designs  (1) 
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Figure  2:  Maximum  initial  plate  deflections 


Longitudinal  bending  of  a  ship  hull  gives  rise 
to  compressive  loads  on  the  cellular  components  of 
a  double  hull.  Hence  the  compressive  Stability  of 
these  components,  in  conjunction  with  lateral 
hydrostatic  presure  loads,  needs  to  be  considered 
in  design.  For  practical  ranges  of  plate  slenderness 
ratios,  local  plate  instability  will  occur  first,  which 
in  turn  affects  the  stiffness  and  strength  of  the 
cellular  beam-columns  spanning  between  transverse 
bulkheads,  in  a  coupled  local-flexural  type  of 
overall  member  failure.  It  follows,  therefore,  that 
the  behavior  of  the  component  plates  must  be 
determined  before  the  overall  behavior  of  the  box 
beam-columns  can  be  studied.  Among  other 
factors,  local  plate  behavior  is  influenced  by  plate 
initial  imperfections,  distortion  from  the  welding 
process,  and  welding  residual  stresses 

This  paper  considers  the  effects  of  these  two 
factors  on  local  plate  strength  and  stiffness.  The 
measurement  of  these  imperfections  is  discussed, 
and  results  from  the  finite  element  analysis  using 
these  measured  imperfections  are  presented. 

INITIAL  PLATE  DEFLECTIONS 

Initial  plate  deflections  are  a  result  of  various 
factors,  such  as  the  transverse  shrinkage  of 
longitudinal  edge  welds,  and  handling  during  the 


fabrication  of  the  welded  structure.  Real  plates  are 
therefore  not  perfectly  flat.  These  kinds  of 
deviations  from  the  flat  surface  (or  residual 
stresses,  see  following  section)  suggest  a  load 
deflection  rather  than  a  buckling  problem.  In 
addition  to  the  amplitude  of  the  initial  deflection, 
the  longitudinal  profile  in  relation  to  the  plate 
aspect  ratio  has  an  effect  on  the  behavior  of  the 
plate.  Carlsen  and  Czujko  (4)  have  shown  that  the 
presence  of  non-elastic  buckling  modes  in  the  total 
modal  content  of  the  imperfections  will  always 
have  a  strengthening  effect  on  a  rectangular  plate. 
However,  for  design  purposes  some  degree  of 
preferred,  elastic  buckling  mode  pattern  in  the 
initial  deflections  must  be  assumed. 

Therefore,  in  order  to  assess  the  effect  of 
initial  deflections  on  plate  strength,  the  initial 
deflections  of  each  component  plate  of  the 
multicellular  box  specimens  were  measured.  The 
length  of  the  specimens  was  182.9  cm  (6')  and  the 
cross  section  consisted  of  single,  double  or  triple 
box  cells.  Square  cells  were  either  457  mm  by 
457  mm  (18"  by  18")  or  914  mm  by  914  mm  (36" 
by  36"),  while  rectangular  cells  were  914  mm  by 
762  ran  (36"  by  30")  or  610  ran  by  762  ran  (24" 
by  30").  The  nominal  plate  thickness  was  9.5  mm 
(3/8"),  giving  rise  to  plate  width  to  thickness  ratios 
of  48,  64,  80  and  96.  Plate  initial  deflections  of 
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Figure  3:  Typical  measured  out-of-flatness  profile  for  a  plate 


the  test  specimens  were  measured  using  an 
aluminum  frame  with  magnetic  pads  at  each  end. 
The  function  of  the  pads  was  to  secure  the  frame 
to  the  longitudinal  edges  of  the  specimen  while  the 
readings  were  taken.  Dial  gauges  were  attached  in 
between  the  pads  and  at  regular  spacings  of  228.6 
mm  or  114.3  mm  (9"  and  4  1/2"),  in  order  to  take 
readings  of  the  plate  surface  profile.  The  reference 
readings  for  the  dial  gauges  were  first  obtained  by 
setting  the  frame  against  a  flat,  milled  surface. 
The  measurements  commenced  at  one  end  of  the 
specimen  at  the  first  station,  where  a  series  of  dial 
gauge  readings  were  taken  across  the  specimen 
width.  Then  the  frame  was  moved  to  the  next 
station,  secured,  and  readings  were  taken  again. 
The  whole  process  was  repeated  until  the  last 
station  at  the  other  end  of  the  specimen  was 
reached. 

The  results  of  the  measurements  are  shown  in 
Figures  2  and  3.  In  Figure  2,  the  maximum 
measured  initial  deflections  as  a  function  of  the 
plate  slenderness  ratio  were  plotted  and  compared 
to  several  equations  for  estimating  maximum  plate 
imperfection  for  a  given  plate,  as  proposed  by 
Carlsen  and  Czujko,  (4)  Faulkner,  (5)  and 
Antoniou  (6).  The  actual  measured  plate 
imperfections  ranged  from  3  to  12  mm  (O.  12  to 
0.47”),  and  were  below  the  predictions  of  the  three 
equations.  Faulkner’s  equation  is  quadratic  and  the 
impetiections  increase  rapidly  with  the  slenderness 
ratio.  Carlsen  and  Czujkos’,  and  Antoniou 's  linear 


equations  show  better  agreement  with  the  measured 
values.  Measured  values  were  also  within  the  U.S. 
Navy  guidelines  for  maximun  allowable  unfairness 
in  ship  plating  as  stipulated  in  MIL-STD-1689SH  - 
Welding,  Fabrication  and  Inspection  of  Ship 
structures. 

Figure  3  shows  the  random  surface  of  one  of 
the  component  plates  of  the  specimens.  Measured 
profiles  were  fitted  to  a  double  sine  series, 
Equation  (1)  below,  using  least  squares  regression 
analysis. 


-  - 


j?y 

b 


(1) 


In  general  the  surface  profiles  are  quite  random 
in  nature,  although  in  a  few  cases,  the  first 
longitudinal  mode  of  deflection,  i.e.,  a  single 
longitudinal  half  wave,  predominates.  For  design 
purposes  it  would  be  impractical  to  measure  actual 
imperfections,  so  some  conservative  procedure, 
such  as  the  proposed  equations  mentioned  above, 
would  be  used. 

WELDING  RESIDUAL  STRESSES 

Residual  stresses  due  to  welding  were 
measured  in  a  single  cell  specimen,  The  stresses 
were  measured  using  the  sectioning  method,  in 
conjunction  with  a  25.4  cm  (10”)  Whittemore 
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_J^<2.5  cm) 


Figure  4:  Residual  stress  test  sample  segment 

gauge;  the  sectioning  method  was  also  used  in 
conjunction  with  strain  gauges.  A  test  segment  of 
44.5  cm  (1’  6")  in  length  was  removed  from  the 
single  cell  specimen  (Figure  4)  and  used  to 
measure  the  residual  stresses  by  cutting  it  into 
strips.  Fillet  welds  were  placed  on  both  sides  of 
the  plates  to  join  them  (Figure  5).  Figure  5  shows 
the  locations  where  Plates  A  and  B  of  the  segment 
were  cut  into  strips.  Plates  C  and  D  had  similar 
cut  patterns.  In  all  plates,  thinner  strips  were  cut 
near  the  welds  where  a  significant  gradient  in 
residual  stresses  was  anticipated.  The  middle  strips 
of  each  plate  were  wider,  where  a  more  constant 
residual  stress  was  anticipated.  Prior  to  cutting  the 
segment  into  strips,  a  reference  strain  measurement 
£q  was  taken.  Strain  gauges  allowed  a  direct  strain 
reading  from  each  gauge.  The  Whittemore  gauge, 
a  mechanical  device  having  a  25.4  cm  gauge 
length,  required  two  small  holes  to  be  drilled  in 
each  strip  in  order  to  take  strain  readings.  These 


Figure  5:  Dissection  marks  for  strips  used  for 
residual  stress  measurements 


holes  were  drilled  at  the  gauge  length  of  25.4  cm, 
centered  in  the  middle  of  the  strip.  After  placing 
the  instrument  and  seating  it  in  a  pair  of  drilled 
holes  in  a  strip,  a  reference  reading  was  obtained 
for  the  strip  by  averaging  three  Whittemore 
measurements.  A  temperature  reference  reading 
was  also  taken,  in  order  to  minimize  the  effects  of 
any  temperature  fluctuation  of  the  steel  during  the 
measuring  period.  A  liquid  coolant  was  applied  to 
the  steel  to  dissipate  any  heat  created  from  the 
cutting  process;  this  way,  any  secondary  residual 
stresses  due  to  thermal  expansion  were  reduced. 
Additional  strain  measurements  £f  were  taken  after 
the  strips  were  cut  The  difference  between  these 
and  the  reference  readings  £q  provided  the  residual 
stresses 

°r«  =  E(*f  ~  £o)  (2) 

where  E  is  Young’s  Modulus,  203.4  GPa  (29,500 
ksi). 

Residual  stress  measurements  were  taken  on 
both  sides  of  all  four  plates.  Figure  6  shows  the 
residual  stresses  measured  on  the  outside  of  the 
plates  using  the  Whittemore  gauge.  The 
Whittemore  gauge  results  show  a  nearly  constant 
compressive  residual  stress  distribution  over  the 
middle  portion  of  each  plate,  with  a  high  value  of 
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+38.3  ksi 
(264.1  MPa) 


Figure  6:  External  Surface  Longitudinal  Residual  Stress  Distribution-Whittemore  Gauge 


PLATE  C 


Figure  7:  External  Surface  Longitudinal  Residual  Stress  Distribution-Strain  Gauges 


tensile  residual  stress  existing  at  the  edges  of  each 
plate  near  the  locations  of  the  welds.  The  residual 
stress  measurements  based  on  the  strain  gauge 
reading  are  shown  in  Figure  7.  Only  a  portion  of 
plates  A  and  D  were  strain  gauged,  as  indicated  by 
the  location  of  the  results  shown  in  this  figure. 


The  strain  gauge  results  agree  closely  with  the 
Whittemore  results,  which  were  taken  at  similar 
locations 

From  the  measured  residual  stresses  of  each 
plate,  an  idealized  rectangular  pattern  of  uniform 
tensile  stresses  up  to  yield  strength  ay  along  the 
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Tensile  Specimen  L2 


Plate 

Average  are 

MPa 

U 

A 

48.3 

3.46 

B 

52.5 

3.74 

C 

64.7 

4.53 

D 

49.5 

3.54 

Table  I:  Values  for  cre  and  rj  from  the 
longitudinal  residual  stress  distribution 


(a) 


(b) 


U _ Jl 

I  Y  33.45'  (84.96  cn)  I  Vi 

'  1.275'  '  1.275' 

(3.24  cn)  (3.24  cn) 

Figure  8:  Measured  and  idealized  pattern  of 
longitudinal  residual  stress,  Plate  A 

longitudinal  edges  of  the  plate  was  calculated, 
balanced  by  uniform  compressive  stress  ore  over 
the  middle  strip.  For  each  plate,  the  compressive 
stress  was  averaged  over  the  middle  strip  of  the 


Figure  9:  Representative  stress  strain  curve  for 
HSLA-80  tensile  coupons 

plate.  The  width  of  the  tensile  edge  strips  was 
then  determined  from  equilibrium  of  the  residual 
forces  (Figure  8).  The  computed  values  of  the 
tensile  width  parameter,  t|,  was  found  to  vary  from 
3.46  to  4.53  for  the  four  plates  (A,  B,  C,  D)  which 
formed  the  cell  of  the  residual  stress  test  segment 
(Table  I).  Faulkner  (5)  has  recommended  values 
of  T|  from  3.0  to  4.5  for  ship  plating. 

FINITE  ELEMENT  ANALYSIS 

The  component  plates  of  the  cellular  box 
specimens  were  analyzed  by  the  finite  element 
(FE)  method.  The  actual  measured  initial  plate 
deflections,  and  idealized  residual  stresses  (where 
r(=3.81),  were  used  in  the  analysis.  The  plates 
were  modeled  with  midsurface  thin  shell  elements 
using  a  large  displacement  formulation  with 
material  plasticity.  Since  the  actual  initial 
geometric  imperfections  were  used,  there  was  no 
plate  symmetry  and  hence  all  the  plates  were 
modeled  whole.  A  bilinear  kinematic  hardening 
material  model  was  used  together  with  the  von 
Mises  yield  criterion.  Based  upon  the  results  of 
tensile  coupon  tests,  the  average  yield  strength  ay 
was  equal  to  600  MPa  (87  ksi)  and  the  strain 
hardening  modulus  Est  was  equal  to  1195  MPa 
(175  ksi).  Figure  9  shows  a  typical  stress  strain 
curve  from  a  coupon  test  performed  on  HSLA-80 
steel. 
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EFF.  OF  IMPERFECTIONS  ON  PLATE  STRENGTH 


Figure  10:  Comparison  of  load  shortening  curves,  b/t  =  45.7 


b 

(mm) 

b. 

t 

2s 

ay 

(n=4) 

b 

Maximum  Strength 

Square 

Plate 

’Perfect’ 

Plate 

Actual 

Plate 

457 

45.7 

0.200 

0.0072 

0.6332 

0.7479 

0.6174 

914 

92.5 

0.091 

0.0072 

0.3946 

0.5548 

0.3825 

Table  11  Comparison  of  maximum  plate 
strengths 

The  residual  stresses  were  entered  directly  into 
each  element  as  initial  residual  strains.  The 
idealized  rectangular  pattern  was  assumed  to 
remain  constant  over  the  length  of  the  plate  except 
near  the  ends,  where  it  linearly  decayed  to  zero 
over  a  length  of  one  plate  width.  Boundary 
conditions  were  taken  to  be  simply  supported  and 
unrestrained  along  the  longitudinal  unloaded  edges, 
and  simply  supported  along  the  loaded  edges. 
Uniform  edge  compressive  loading  was  applied  by 
means  of  prescribed  uniform  shortenings  of  the 
plate  from  one  end,  and  the  total  load  was  obtained 


from  the  reactions  at  the  other  end. 

Figures  10  through  13  show  the  results  of  the 
finite  element  analyses  for  two  rectangular  plates 
with  nominal  slenderness  ratios  b/t  equal  to  48  and 
96,  respectively.  The  nominal  aspect  ratios  a/b  are 
equal  to  4  and  2.  Figures  10  and  11  compare  the 
load  shortening  curves  for  the  actual  plate,  the 
corresponding  "perfect"  plate  with  no  residual 
stresses  and  initial  deflections,  and  also  the 
corresponding  square  plate  with  the  same  residual 
stresses  and  amplitude  of  initial  deflection  as  the 
actual  plate.  The  imperfection  profile  for  the 
square  plate  was  taken  to  be  single  half  waves  in 
both  directions.  In  both  plates,  it  can  be  seen  that 
plate  imperfections  have  a  marked  influence  on  the 
stiffiess  and  strength  of  the  plates  the  influence  is 
also  greater  for  the  more  slender  plate  compared  to 
the  less  slender  (stockier)  plate.  Table  11  lists  the 
maximum  strengths  of  the  plates,  where  it  is  seen 
that  the  percentage  strength  reduction  due  to  the 
combined  effects  of  initial  deflections  and  welding 
residual  stresses  is  about  30%  for  the  slender  plate, 
compared  to  1690  for  the  stockier  plate.  The 
imperfections  also  reduced  the  gradient  of  the 
unloading  portion  of  the  curve  in  the  case  of  the 
stockier  plate,  as  shown  in  Figure  10.  The  load 
shortening  behavior  of  rectangular  plates  is  often 
conservatively  assumed  to  be  similar  to  that  of  an 
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EFF.  OF  IMPERFECTIONS  ON  PLATE  STRENGTH 


Figure  11:  Comparison  of  load  shortening  curves,  b/t  =  92.5 


Figure  12:  Deflected  shape,  b/t  =  45.7 


Figure  13:  Deflected  shape,  b/t  =  92.5 


equivalent  square  plate.  The  figures  show  that  the 
square  plate  curves  are  less  stiff  than  the  actual 
rectangular  plate  curves;  peak  strengths  also  do  not 
differ  significantly.  Figures  12  and  13  show  the 
deformed  shape  of  the  plates  where  the  number  of 
half  waves  formed  is  clearly  seen  to  be  4  and  2. 

The  maximum  strengths  of  all  the  separately 
analyzed  plates  of  the  test  specimens  (see  Table 
III)  are  plotted  and  compared  with  seeral  design 


equations  in  Figure  14.  The  design  curves  are  the 
Frankland  (7)  curve  which  is  used  by  the  U.S. 
Navy,  and  the  Faulkner  (5)  curve.  The  latter  is 
plotted  for  two  cases;  with  and  without  residual 
stresses.  The  elastic  buckling  curve  is  also  shown 
for  comparison.  Frankland ’s  equation  gave  good 
predictions  at  the  higher  slenderness  ratios,  but 
overpredicted  somewhat  for  the  stockier  plates 
(b/t=48  and  64).  By  contrast,  all  the  data  points 
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STRENGTH  OF  RECTANGULAR  PLATES 


-  AD1NAFEM  plates 


Figure  14:  Comparison  of  FEM  plate  results  with  plate  design  equations 
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2.461 
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4.147 
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2.461 

0.6426 

4.249 

0.4549 

2.553 

0.6334 

4.249 

0.4676 

2.553 

0.6345 

4.927 

0.4050 

2.553 

0.6371 

4.927 

0.4002 

3.426 

0.5283 

4.981 

0.3983 

3.426 

0.5417 

4.981 

0.4000 

3.426 

0.5400 

4.981 

0.4002 

3.371 

0.5298 

5.029  ’ 

0.3984 

3.371 

0.5339 

5.029 

0.4008 

3.371 

0.5322 

5.116 

0.3876 

4.055 

0.4690 

5.116 

0.3906 

4.060 

0.4758 

5.116 

0.3922 

4.114 

0.4711 

Table  III:  Maximum  strengths  of  rectangular 
plates  by  finite  element  analysis 


were  above  the  Faulkner  stress-free  plate  equation. 
The  Faulkner  equation  for  welded  plates  is  seen  to 
give  very  conservative  predictions  compared  to  the 
finite  element  results.  The  results  indicate  that  the 
Frankland  equation  is  suitable,  while  the  Faulkner 
equation  seems  to  be  the  best  for  estimating  plate 
ultimate  strengths. 

CONCLUSIONS 

Initial  imperfections  have  been  shown  to  affect 
the  behavior  of  welded  plates.  The  actual 
measured  initial  deflections  of  the  welded 
component  plates  of  a  prototype  double  hull 
structure  were  in  the  range  of  3  to  12  mm  (0.12  to 
0.47"),  and  within  specifications.  Measured 
residual  stresses  in  the  component  plates  of  the 
welded  module  showed  the  typical  pattern  of 
narrow  edge  strips  of  high  tension,  and  a  wide 
middle  strip  of  low,  fairly  uniform  compressive 
stresses.  The  maximum  strengths  of  welded  plates 
were  reasonably  predicted  by  the  Frankland  and 
Faulkner  stress-free  equations.  Results  also 
showed  that  equivalent  square  plate  load  shortening 
behavior  can  be  used  to  represent  the  behavior  of 
the  corresponding  rectangular  plate. 
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